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Abstract 



We propose a simple test of underground neutrino data to discrimi- 
nate amongst neutrino oscillation models. It uses the asymmetry between 
downward-going events and upward- going events, for electron and muon 
events separately. Because of the symmetry of typical underground detec- 
tors, an asymmetry can be compared with calculations with little need for 
the intermediary of a simulation program. Furthermore, we show that the 
various oscillation scenarios give rise to dramatically differing trajectories of 
asymmetry versus energy for muons and electrons. This permits a clean dis- 
tinction to be drawn between models. 
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I. INTRODUCTION 



The atmospheric neutrino anomaly has been with us now for more than ten years. It 
was discovered serendipitously in the largest detectors built for nucleon decay searches, 1MB 
|l|] and Kamioka 0, and confirmed recently in Soudan [Q. It appeared as a deficit in the 
ratio of muon-like to electron-like neutrino interactions within the fiducial volumes of these 
massive Cherenkov detectors. The atmospheric neutrinos in the range of 1 GeV of energy 
originate high in the atmosphere from the decay of pions and muons. The ratio of the z/^ 
flux to z/g flux is thus closely constrained by well known decay kinematics to be in the ratio 
of 2:1 (actually, 2.2:1). The experiments have found that the apparent ratio of observed flux 
flavor ratio to expected flux flavor ratio is closer to R = 0.60 ± 0.05. The result is usually 
presented in terms of this ratio of ratios because the absolute neutrino fluxes are not very 
well known (~ 20%, or perhaps even worse), whereas the 0(z/e)/0(z/^) is known to several 
percent. The simplest explanation for the anomaly seems to be the one originally suggested, 
namely neutrino oscillations [Q. 

Until the present time the size of the underground detectors has limited the number 
of events to a few hundred, and the energy of those events to less than about 1.5 GeV. 
Previous data have not shown conclusive zenith angle dependence as expected from neutrino 
oscillations. At energies of a few GeV however, we must begin to observe angular effects, if 
neutrino oscillations are the cause of the anomaly. Indeed, the flrst public presentations of 
data from the new massive Super-Kamiokande experiment do seem to indicate some angular 
variation p. 

Our purpose in this note is to set the stage to interpret future data in terms of discrim- 
inating amongst the many scenarios which have been constructed to explain one or more of 
the three outstanding "problems" in neutrinos: the solar neutrino problem 0, the LSND 
effect 0, and the presently considered atmospheric neutrino anomaly. 

We deflne an asymmetry in direction as simply 



where D are the number of downward-going events and U are upward-going events, for 
each of muon neutrino events and electron neutrino events. We assume the detector to be 
up/down symmetric, and the data set to be free of significant (presumably downward-going) 
background and crossover between muon and electron type of charged-current events. 

We ignore the effect of Uj. in our calculations, as well as the contamination of events by 
neutral current interactions. The Ur charged-current cross section is sufficiently small at the 
energies discussed herein that it makes a negligible contribution (to both muon and electron 
events). Thus for present considerations, oscillations between i/^ and Ur are indistinguishable 
from oscillations between z/^ and a new sterile neutrino species. Neutral currents are a small 
fraction of events, and in any case should show no asymmetry (except in the sterile neutrino 
case) . 

II. CALCULATIONS 

We now calculate this asymmetry for a wide variety of oscillation scenarios. We consider 
the following cases: 

(A) two- flavor mixing (a) z/^ — Uj. and (b) — Ve'-, 

(B) three-flavor mixing with a variety of choices for mass-mixing parameters; 

(C) sterile maximal mixing; and 

(D) massless neutrino mixing. 

It will be shown that these various scenarios lead to very different predictions in the sign, 
magnitude or the energy dependence of the asymmetry and hence are (relatively) easy to 
distinguish from one another once there is sufficient data available. 

A. Two-Flavor Mixing 

a. Ufj, — Ur Mixing. This is the simplest case possible Q]. The z/g flux is unaffected and 
the flux modified as 
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where P^^ = 1 — sin^ 29 sin^ . (2) 

Ah/ 



Sm^L 



Hence A^, is essentially zero and independent of energy. At low enough energies P^ 
1 — I sin^ 29 and is independent of L and hence ^ at low energies. At high energies; 
L/E is negligible for down u'^s and hence 

Ni = at;. (3) 

For upward-going u'^s P^^ ~ (1 — | sin^ 29) and NJ^ ^ (1 — | sin^ 2^) A'° and hence: 

_ l/2sin^2g 
l + l/2sin2 2^ 

has a maximum of 1/3 when sin^ 20 = 1. Note that at high enough energies will 
come back asymptotically to zero. 
b. — z/g Mixing In this case 0] 

N, = N'^iP +r(l-P)) 

Ar, = Ar°(P+(l/r)(l-P)) (5) 

where P = P,,,, = Ppp, as in Equation (2), and r = N^/N^. A gain at low energies Ae = 
A^ ^ 0. At high energies 

Ni = Nl, Ar; = ArJ (p + r(l-P)) 
AT,'^ = AtO, at; = AT" (P + (l/r)(l - P)) 

^ (l-P)-l/r(l-P) 



(l + P) + l/r(l-P) 
For P = 1/2 we get the limiting values: 



1 — r 



and 



Recall that r ~ 0.45 at low energies, decreases to 0.3 at -Ej, ~ 5 GeV and eventually 
becomes negligible. Note that and Ae always have opposite signs in this case. 



B. Three-Flavor Mixing 

There are two ways to account for all three neutrino anomalies (solar, atmospheric and 
LSND) with just three flavors. In one, due to Cardall and Fuller ^j, a single Sm"^ is expected 
to account for both the atmospheric low energy anomaly as well as the LSND observations. 
In particular, 0.3 eV"^ ~ ^^31 ~ ^^32 ^^12 ~ 10~^ eV"^, with large — 1/^- mixing. In 
this case the resulting probabilities are very similar to the two-flavor — Vj. mixing, but 
with a large 5m? of 0.3 eV'^. As a result of the large 5m? very little zenith angle dependence 
or asymmetry (neither A^, nor A^) is expected. 

The other three- flavor solution, due to Acker and Pakvasa accounts for both solar 
and atmospheric anomalies with a single and with large mixing between z/g and i/^, the 
mass pattern being 2 eV'^ ~ ^^li ~ '^'^32 >> '^^12 ~ 5 x 10^'^ eV"^ . The probabilities in 
this case are essentially identical to the two-flavor — Vf, case with large mixing. In both 
of these scenarios, it is possible to have nearly degenerate neutrinos with cosmologically 
significant total mass. 

There are also other scenarios with three-neutrino mixing with a wide range of mixing 



patterns |T0|. In general, we expect them to yield asymmetries which will interpolate between 
the two limiting cases of z/^ — and — z/g mixing. 

We consider as a unique and interesting example, the maximal mixing proposal of Harri- 
son, Perkins and Scott |]1T|. The assumption is that 5m?^2 ~ '^"^31 >> ^^12 ~ 10~^^eV^^ and 
^mg^ is in the range of 10^^ — 10~^ eV"^ . With the assumed maximal mixing, the probabilities 
for L/ E in the appropriate atmospheric range are given by 

P^^. = Pee = l- 8/9 sin^ 



P,e = Pe,{= Per = P,r) = 4/9 sin^ (9) 

The expected asymmetries can be easily written down 

(1 - P,,) - rP,, (1 - Pee) - l/rPe, 

' (l-P,,)+rP,e' ^ (1 + Pee) - 1/rPe, ^ ^ 

Hence both and A^, are small at low energies and at high energies P^^ = P^^ = 1 and 
Pe'j, = 0, whereas PJi^ = P^^ = 5/9 and P^^ = 2/9. And so 

^ _ 4/9 -r 2/9 _ 2 - r 
^ ~ 14/9 + r 2/9 ~ 7 + r 

and A - - - -i^ fin 

^"""^ - 14/9 + l/r(2/9) " 7r + l ^^^^ 



C. Sterile Maximal Mixing 



In the scheme of Foot and Volkas |T2| , z/^ mixing maximally with a new sterile i/^/ with 
a Srn'jj^^, ~ 5 x 10^'^ eV"^ accounts for the low-energy atmospheric neutrino anomaly. Hence, 
the muon asymmetry is identical to the one in the case of t'^ — Ur oscillations, as discussed 
above. In addition, Ug mixes maximally with a sterile z/g', and when Sml^, is in the range 
of lO'^eV"^ electron-neutrinos will also oscillate and get depleted. The resulting electron 
asymmetry is strikingly different from the case of z/^ — oscillations in having the opposite 
sign, which makes it unique and easy to distinguish. Ae and will differ only in slightly 
different energy (or L/E) dependence but be otherwise similar and always have the same 
sign. 



D. Massless Neutrinos 

If neutrinos are massless, there can still be mixing and oscillations. Two possibilities have 
been considered in the literature. One is the case where different flavors couple differently 
to gravity and the other is a breakdown of Lorentz invariance where each particle may 
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have its own maximum speed [|^. The oscillation phenomenology is identical for both cases. 



The survival probability in the two-flavor limit is given by 

P^f, = Pee = l-sm^29sm^ (Isv El) (12) 



2 

where 6 is the flavor mixing angle, and 6v is the small parameter characteristic of violation 
of equivalence principle or Lorentz invariance. Note the strikingly different dependence 
on L and E: L x E instead of L/E. Remarkably, an allowed choice of parameters is 



able to account for both solar and atmospheric neutrinos [|T5l: sin 26^ ~ 0.8 to 1 and 
- ■ 10~^ — 10~^(fcm — GeV)~^. The expressions are the same as in the — Ue case except 



that sin^ ^ is replaced by sm'^{6v/2LE). As a result, the roles of low and high energy are 
reversed. The asymmetries and Af, become rather small at high energies; at low energies 
they are given by Eq.(12). 



III. NUMERICAL RESULTS 

We have performed numerical calculations of the models discussed above. They are 
explicitly: 

1. Simple two-flavor oscillations between and [Q. The example is for Sm^ = 
0.005 eV^, sm^29 = 1. 

2. Two-flavor oscillations between i/^ and z/e with the same parameters as above. (The 
Acker-Pakvasa ^ scheme leads to the same result). 

3. Three-flavor mixing a' la Cardall- Fuller M. 



4. Three-flavor maximal mixing scheme of Harrison-Perkins-Scott [O . 



5. Sterile maximal mixing of Foot-Volkas W2 



6. Massless neutrino mixing, where we take 6v/2 ~ 10 ^ {km — GeV) ^ [ITB 
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In Figure |l| we show results for and as functions of energy from more detailed 
calculations. We calculated energy spectra between 0.2 and 5.0 GeV for a detector with 
an exposure of 22 kiloton-years (approximately one year of Super-Kamiokande data). We 
use the Bartol flux model, and a simple quark model for the charged-current cross section, 
and assume a perfect detector ||TB||. Detailed calculations for a particular instrument will 
of course vary, but the asymmetry will change little, the general behavior illustrated being 
insensitive to the details. 

We show the trajectories of versus Ag in Figure ^ for the six models. Note that there 
are small asymmetries at low energies, due to the inhomogeneity of the earth's magnetic 
field, as incorporated into the atmospheric flux model we employ |jl6|. It is clear that with 
good statistics all scenarios can be clearly distinguished by both energy dependence and 
relative signs of and Af,. In particular it is noteworthy that the flrst model, currently 
seemingly favored in preliminary reports from Super-Kamiokande 0, stands out distinctly 
from all others. It is straightforward to plot the expected asymmetries in other scenarios or 
different choices of parameters. 



IV. CONCLUSIONS 

In the foregoing we have presented a case for employing the up-to-down asymmetry of 
neutrino interactions in underground detectors as a discriminator for some of the many 
neutrino oscillation schemes which have been discussed as solutions to various combinations 
of the current three neutrino puzzles (solar, atmospheric and LSND). The asymmetry has the 
virtue that it can be calculated directly from data (using only particle identiflcation, energy 
and direction), without aid of simulation programs. It is self-normalizing and independent 
of flux model calculations, and tests electron and muon data separately. 
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FIGURES 




Model 1 





Model 2 





Model 3 




Model 4 



Model 6 



Model 5 

FIG. 1. The muon (solid lines) and electron (dashed lines) asymmetries versus energy (GeV), 
for the 6 oscillations models considered herein (see text). 
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FIG. 2. The trajectories of the muon asymmetry and electron asymmetry for the 6 oscillation 
models considered herein. The arrowheads point in the direction of increasing charged lepton 
energy, which ranges from 0.2 to 5.0 GeV in these calculations. 
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